Localized waves whithout the existence of extended waves: 
oscillatory convection of binary mixtures with strong Soret effect 
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Spatially confined solutions of traveling convection rolls are determined numerically for binary 
mixtures like ethanol- water. The appropriate field equations are solved in a vertical crossection of the 
rolls perpendicular to their axes subject to realistic horizontal boundary conditions. The localized 
convective states are stably and robustly sustained by strongly nonlinear mixing and complex fiow- 
induced concentration redistribution. We elucidate how this enables their existence for strongly 
negative separation ratios at small subcritical heating rates below the saddle-node of extended 
traveling convection rolls where the quiescent fiuid is strongly stable. 
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Many nonlinear dissipative systems that are driven suf- 
ficiently far away from thermal equilibrium show selfor- 
ganization out of an unstructured state: A structured 
one can appear that is characterized by a spatially ex- 
tended pattern which retains some of the symmetries of 
the system [1]. Some systems form also spatially con- 
fined structures [2] . In particular localized traveling wave 
(LTW) structures that are embedded at subcritical driv- 
ing rates in the stable featureless surrounding and that 
compete with the subcritically bifurcated extended trav- 
eling wave (TW) pattern have been explored in experi- 
ments [3], numerical simulations of the appropiate field 
equations [4, 5], and via various model approaches [6, 7]. 
The coexistence of the TW state and the quiescent state 
is a prerequisite for interpreting LTWs as spatial connec- 
tions of the former to the latter via two fronts. However, 
this picture coming from simple cubic-quintic amplitude 
equations is too simple for convection in binary mixtures: 
there the strongly nonlinear TWs cannot be described 
by a power-law expansion in one amplitude [8]. Here we 
show that LTWs are - independent of TWs - separate 
nonlinear states. 

We predict and describe in quantitative detail stable 
subcritical LTWs in mixtures with strong Soret effect at 
small heating rates where no extended TW whatsoever - 
stable or unstable - exists [9]. Furthermore, we elucidate 
that the concentration field and its current are the key to 
understand how LTW convection is sustained when the 
quiescent conductive state is strongly stable. The latter 
is solutally stabilized since the Soret coupling between 
deviations ST of temperature and SC of concentration 
from their means increases (decreases) the lighter com- 
ponent in cold (warm) fluid regions thereby reducing the 
buoyancy. The most important effect is a flow-induced 
concentration redistribution within the LTW: It is the 
strongly nonlinear mixing by the LTW flow with large 
Peclet numbers that reduces the Soret-generated, stabi- 
lizing concentration gradients and that increases locally 
the driving buoyancy force to levels which suffice to cause 
convection in well mixed fluid with strongly reduced con- 



centration variations [12]. 

We have solved the field equations [13] for convection 
in binary mixtures like ethanol-water for Lewis number 
L=0.01, Prandtl number a— 10, and separation ratios [1] 
-0.25 <-(/'< -0.4 with a finite-difference method [14] 
in a vertical cross section through the convection rolls 
perpendicular to their axes. To simulate ID patterns 
arising in narrow annular channels we applied laterally 
periodic boundary conditions at a; = 0, F with F up 
to 160 times the layer height to accomodate the largest 
LTWs. Spatially extended TW states of different wave- 
length X — 2TT/k were obtained with F = A. Horizontal 
boundaries at top and bottom, z = ±1/2, are realistic, 
i.e., no slip, perfectly heat conducting, and impermeable. 
Our control parameter r — R/R^ measuring the thermal 
driving is the Rayleigh number R reduced by the critical 
one i?°=1707.762 for onset of convection in a pure fluid. 
To measure the lateral variation of the mixing in oscilla- 
tory convective flow we use the mixing number M{x) 
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Here overbars imply vertical averaging. Brackets denote 
temporal averaging, (/) = {f{x -t- w^i, z, t)), over one os- 
cillation period in the frame comoving with the slow drift 
velocity of the LTW envelope. In a perfectly mixed 
fluid M vanishes while M — > 1 in the conductive state 
[denoted by the subscript cond in Eq.(O.l)] with its large 
Soret-induced vertical concentration gradient. 

In Fig. 1 we show how an increasing Soret cou- 
pling strength changes the bifurcation properties of LTW 
width I, maximal vertical flow velocity Wmax^ frequency 
Lo, and LTW drift velocity Vd- The driving interval 
{fmim Tmax) with Stable Stationary LTWs moves upwards 
and its extension increases with \il>\. At its lower end 
LTWs are for all ^ narrow pulses of universal structure 
and I of about 5. At the upper end I seems to diverge. 
LTWs with, say, I < 7 behave differently than the broad- 
width LTWs (cf. Fig. 2) with an extended TW like center 
part between the two limiting fronts. The change in the 
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LTW bifurcation diagrams of lo and versus r in Fig. 1 
reflects this difference. Our simulations indicate that the 
LTWs of Fig. f are uniquely selected, monostable con- 
fined solutions [15] and robust against small perturba- 
tions. 

With increasing |?/;| TW saddle-node location and 
oscillatory threshold Tosc move much faster to larger r 
than the lower LTW band edge Vmin- At ^=-0.08 (not 
shown here) one has < rose < I'min so that there 

the narrow LTW pulses that have been investigated ex- 
tensively lie in the linearly convectively unstable param- 
eter region (see, e.g., [13] for discussion and references). 
However, already at ^'=-0.4 almost all LTWs appear be- 
low the TW saddle-node, i.e., ahead of the TW nose of 
Fig. 3. Thus, stable LTW convection is driven here at 
heating rates r for which any extended TW convection is 
impossible. Even at ■!/'=-0.25 the relevant LTW parame- 
ter combinations of u), r, kpiateau (big bullets in Fig. 3) lie 
outside and ahead of the TW bifurcation surface. For all 
■0 investigated here the LTW width seems to diverge with 
increasing r when its frequency and plateau wave number 
kpiateau (big bullets in Fig. 3) approaches a,n ui — r — k 
combination of a TW. 

LTWs are sustained by the following complex, large- 
scale concentration redistribution process. While trav- 
eling from tail to head LTW rolls increase their lateral 
concentration contrast, i.e., M and with it Vp, and A [16] 
(cf. Fig. 2). Positive "blue" (negative "red") concen- 
tration deviation from the global mean is sucked from 
the top (bottom) boundary layer into right (left) turn- 
ing rolls as soon as they become nonlinear under the 
trailing front. This happens when the vertival velocity 
w roughly exceeds Vp [left arrow in Fig. 2(b)] so that 
regions with closed streamlines appear [8, 13]. Within 
them "blue" ("red") concentration is transported pre- 
dominatly in the upper (lower) part of the layer to the 
right. Mean concentration, on the other hand, migrates 
mostly to the left along open streamlines that meander 
between the closed roll regions and that follow the bor- 
der line between green and yellow in Fig. 2(a). The time 
averaged current of 5C [green lines in Fig. 2 (e)] reflects 
the mean properties of this transport. Since "blue" and 
"red" ("yellow") is transported away from (towards) the 
left trailing front mean concentration accumulates there 
and causes a strong drop of M(x). By the same token 
the leading front's concentration variations and with it 



M (x) are strongly increased even beyond the conductive 
state's values. Thus, unlike TWs LTWs do not reach a 
balance between SC injection and advective mixing and 
diffusive homogeneization on a constant level of small M. 
Rather LTW rolls collaps under the leading front when 
Vp has grown up to w [right arrow in Fig. 2(b)]. There- 
after concentration is discharged and sustains a barrier 
of (SC) ahead of the leading front. When the mixing in- 
creases with Wfnax said r so does I since the slower growth 
of 6C and M to the higher levels necessary for the lead- 
ing front's transition to conduction requires longer and 
longer propagation lengths for the rolls. 

The lateral redistribution causes the mean convectively 
generated C-profile [green line in Fig. 2(d)] to extend 
significantly further into the conductive region than the 
mean convective temperature field (red line). Thus, the 
buoyancy (b) [black line in Fig. 2(d)] is determined in the 
front regions predominantly by the concentration field. 
This explains (i) the decrease of buoyancy below conduc- 
tion levels ahead of the leading front with the associated 
restabilization of conduction there and (ii) the increase 
of (6) out of the conductive state under the trailing front 
and its strong overshoot [cf. Fig. 2(d)] over the bulk 
enabling convection growth. 

That the driving motor for LTW convection is located 
under the trailing front where (b), M,Vp, X are extremal 
can also be seen from the relaxation behavior after chang- 
ing r: Not only {b),M, and Vp but also ui and Vd appro- 
priate for the new r are almost instantaneously realized 
locally under the trailing front while the relaxation of 
the leading front involving also diffusive processes takes 
much longer time, cf.. Fig. 4. Thus there lu and Vd is 
largely selected at the trailing front while the large-scale 
lateral concentration transport and redistribution deter- 
mines the leading front's properties and its location, i.e., 
the LTW width. 

To summarize: We have found, analyzed, and ex- 
plained uniquely selected stable subcritical LTWs in mix- 
tures with strong Soret effect at small heating rates where 
extended TWs are not possible, i.e., below the saddle 
node bifurcation of the latter. These strongly nonlinear 
LTWs are self-consistently sustained by a concentration 
redistribution such that flow-induced mixing increases 
locally the buoyancy to levels that suffice to drive well 
mixed fluid flow. 
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FIG. 1: Bifurcation properties of LTWs (symbols) and TWs (lines) for different separation ratios ip: (a) Pull width I of LTWs 
at half maximum of the envelope of the vertical velocity field w [blue line in Fig. 2(b)]. For small filled symbols I kept on 
growing slowly, (b) Maximal vertical flow velocity Wmax- (c) Frequency; for LTWs in the frame comoving with Vd- (d) Drift 
velocity Vd of LTWs. Phase velocity is always to the right, positive, and much larger. Thick lines in (b, c) denote TWs with 
saddle-node wave number kj^ ~ tt. Thin lines are TWs with plateau wave number kpiate.au of the last LTW before the I oo 
transition from LTW to TW. The frequency of the latter is determined in the frame comoving with the last LTW. Unstable 
TWs (dashed lines; determined with a control method) bifurcate subcritically with large Hopf frequency u>H{k) at rosc{k) out 
of the conductive state and become stable (solid lines) at the saddle-node rj^ (k) when lateral periodicity is imposed with 
r = A = 2n/k. 



FIG. 2: Broad LTW of 1=17.4: (a) Snapshot of concentration deviation 5C from global mean (light green/yellow) in a vertical 
cross section of the layer, (b) Snapshots of lateral wave profiles at midheight, z=0, of SC (green), vertical velocity w (blue), 
and its envelope. At the arrows Wmax = Vp. (b) Mixing number M (green), Eq. (0.1) and phase velocity Vp (black) of nodes 
of w{z = 0) in frame comoving with Vd- The variation of A(a::) = 2n Vp{x)/uj is the same since the LTW frequency w is a 
global constant. Thin (thick) bullet marks smallest (plateau) wavelength for discussion, (d) Time averaged deviations from 
the conductive state at z=-0.25 for concentration (green), temperature (red), and their sum ({&)) measuring the convective 
contribution to the buoyancy, (e) Streamlines of time averaged concentration current (J) = {u6C — L^{6C — tpST)) (green) 
and velocity field (u) (blue). The latter results from (b) and documents roll shaped contributions of (u)(5C) to (J) under the 
fronts and the associated (SC) redistribution. Thick blue and green arrows indicate (u) and transport of positive dC (alcohol 
surplus), respectively. Thus, in the lower half of the layer negative 6C (water surplus) is tranported to the right. Parameters 
are V'=-0.35, r=1.346. 



FIG. 3: TW and LTW oscillation frequencies over the k — r plane. Extended TWs (nose shaped surface) start as unstable 
solutions with Lonik) at upper thick line - its projection onto the k — r plane is the threshold curve rosc{k). Dent in upper 
nose part is related to appearence of closed w-streamlines when Wmax — Vp causing (5C-anharmonicities and the breakdown of 
small-amplitude expansions [8, 13]. With decreasing ui TWs become stable in periodic systems of F = 27r/fc at the saddle-node 
line (dotted outermost perimeter of the nose) located at {k). For large F TWs become stable at the slightly larger Eckhaus 
boundary r'j^^ [k) [18]. TWs with kj^ ~ tt of the tip of the nose (plateau wave number kpiateau of the last LTW before its 
I —f oo transition) are marked - as in Fig. 1 - by thick (thin) curves on the nose surface. Thin dotted horizontal lines denote 
LTWs: they contain for fixed u) a broad band of wave numbers extending from the largest ones, kmax — 4.3 (small bullets as in 
Fig. 2), to small ones outside the plot range. Plateau wave numbers kpiateau are marked as in Fig. 2 by big bullets that appear 
only where / is sufHciently large. 



FIG. 4: Thick full (dashed) lines show front relaxation over a limited time interval after changing the driving from r=1.41 to 
1.43 (1.38). Thin gray (black) lines are world lines of w-nodes (front positions) of the original LTW (r=1.41, ip=-OA, 1= 20). 



